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Abstract: We argue that the recently discovered integrabihty in the large- 
CFT/AdS system is equivalent to diffractionless scattering of the corresponding hid- 
den elementary excitations. This suggests that, perhaps, the key tool for finding 
the spectrum of this system is neither the gauge theory's dilatation operator nor the 
string sigma model's quantum Hamiltonian, but instead the respective factorized S- 
matrix. To illustrate the idea, we focus on the closed fermionic 5u(l|l) sector of the 
A/" = 4 gauge theory. We introduce a new technique, the perturbative asymptotic 
Bethe ansatz, and use it to extract this sector's three-loop S-matrix from Beisert's 
involved algebraic work on the three-loop su(2|3) sector. We then show that the 
current knowledge about semiclassical and near-plane-wave quantum strings in the 
su(2), su(l|l) and s[(2) sectors of AdS^ x is fully consistent with the existence of a 
factorized S-matrix. Analyzing the available information, we find an intriguing rela- 
tion between the three associated S-matrices. Assuming that the relation also holds 
in gauge theory, we derive the three-loop S-matrix of the s[(2) sector even though 
this sector's dilatation operator is not yet known beyond one loop. The resulting 
Bethe ansatz reproduces the three-loop anomalous dimensions of twist-two operators 
recently conjectured by Kotikov, Lipatov, Onishchenko and Velizhanin, whose work 
is based on a highly complex QCD computation of Moch, Vermaseren and Vogt. 
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1. Introduction and Conclusions 

There is mounting evidence that four-dimensional strictly planar A/" = 4 Yang-Mills 
gauge theory is integrable. Likewise, there are strong indications that free IIB su- 
perstring theory on the curved space AdS^ x is also integrable. This means that 
there is hope that the spectrum of both theories might be exactly computable. If 
true, the AdS/CFT duality conjecture becomes falsifiable: Either the spectrum of 
the two models agrees, or it does not. This is good news, since, according to scientific 
tradition, falsifiability is an important feature of a theory, to be distinguished from 
creative speculation. 

The first, crucial hint that conformal Af = 4 gauge theory might be integrable 
in the planar limit was discovered in a beautiful paper by Minahan and Zarembo 
[0]. There it was shown that the spectrum of A/" = 4 conformal operators may be 
obtained by diagonalizing an integrable quantum spin chain. The observation was 
initially restricted to the subset of scalar operators at one loop. However, shortly 
after, strong evidence was found that integrability extends to higher loops 0, and. 
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at least at one loop, to the full set of A/" = 4 operators 0. Integrable structures 
have appeared before in planar QCD, starting with the pioneering work of Lipatov 
[^] (for a comprehensive review, see P). There it was always considered to be due 
to approximate, hidden symmetries, which can lead to useful information on high 
energy scattering and QCD anomalous dimensions. In A/" = 4 0, Q the goal is 
more ambitious. One would like to exactly solve a four-dimensional gauge theory. 

In a parallel, intertwined development integrable structures were also observed 
in the worldsheet theory of strings on AdS^ x S^. First, there were some hints that 
the relevant classical coset sigma model is integrable ||^. Later it was shown in detail 
how finite-dimensional reductions of the classical sigma model lead to integrable 
classical mechanics systems of Neumann type The reduction may be undone 
by the Backlund transformation, and the infinite tower of commuting charges of 
the original sigma model is recovered [§. Subsequently, in an important paper by 
Kazakov, Marshakov, Minahan and Zarembo, it was demonstrated how the spectrum 
of finite-gap solutions of the classical sigma model may be obtained by the classical 
inverse scattering method [Q. In [IC] first steps toward finding the quantum spectrum 



of the sigma model were taken, and an approximate S-matrix was proposed. Very 
recently, the Hamiltonian of the classical bosonic string propagating on AdS^ x 5*^ 
was shown to be integrable by constructing, in a special gauge, the corresponding 
Lax representation |Tl|. 

How should one proceed toward, for one, exactly solving the planar "CFT", 
namely the large limit of A/" = 4 gauge theory, and, secondly, the worldsheet 
theory of strings on AdS^ x Let us assume, as a working hypothesis, that both 
models are indeed not only approximately, but in fact completely integrable. If true, 
proceeding with the construction of the gauge theory's dilatation operator, loop- 
by-loop, or meticulously quantizing the full sigma model, including fermions, might 
prove too hard. In fact, integrable systems do not always have simple Hamiltonians. 
To the contrary, the Hamiltonian might be quite intricate in order to realize the 
subtle hidden symmetries responsible for integrability. 

Is there a simpler object which encodes the spectrum of an integrable quantum 
system? In fact there is: The S-Matrix. Quantum integrability is deeply tied to 
the concept of diffractionless, factorized scattering^. It means that the elementary 
excitations of a quantum many-body system interact only through a sequence of 
two-body scattering processes which may lead to the exchange of quantum numbers 
and momenta, but do not alter the magnitudes of the latter. This is the next-best 



^We are discussing in this paper an internal S-Matrix describing the scattering of elementary 
excitations on a lattice hidden inside the trace of of gauge invariant composite local operators. It 
should not be confused with the external S-matrix of A/" = 4 which refers to multi-gluon amplitudes 
in four-dimensional space time. Recently dramatic progress was also achieved in this direction, see 
p^ , and references therein and thereafter. It would be exciting if a relation between the internal 
and the external S-matrix could be found. 
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thing to a free system! Interaction does take place, but scattering only leads to a 
permutation of a fixed set of momenta (and, possibly, of quantum numbers) . 

Here we would like to propose that one should try to directly construct the S- 
matrix of integrable gauge and string theory. Hopefully, the two will then agree. 
Unfortunately we do not have a very concrete proposal about how to go about it in a 
direct fashion. However, since the AdS/CFT system possesses a huge amount of sym- 
metry, and, in particular, apparently an infinite number of hidden conserved charges, 
we may certainly hope that methods to "bootstrap" the S-matrix will eventually be 
found. 

As a more modest first step, we will analyze part of what is known about the 
gauge theory's dilatation operator, and about the spectrum of AdS strings, in order 
to get a glimpse at how the S-matrix might look like. We think that the results are 
encouraging. For simplicity we will restrict ourselves, in this paper, to three simple 
two-component sectors of the full superconformal, psu(2,2|4) symmetric, system. 
These are a compact su(2) subsector (two bosons), a compact su(l|l) subsector (one 
boson, one fermion), and a non-compact sl{2) subsector (two bosons). 

Let us recall what is currently known about the A/" = 4 gauge theory's dilatation 
operator and its Bethe ansatz. The latter may be constructed once the S-matrix 
is known. At one loop, the dilatation operator was derived for the complete set of 
psu(2,2|4) operators |13[- It was shown to be integrable in 0, and the associated 
Bethe ansatz was constructed. Beyond one loop, we have Beisert's two- and three 
loop dilatation operator in the "maximally compact" subsector 5u(2|3) [0. This 
includes the su(2) (here it was initially found in 0) and su(l|l) sectors, but not 
the s[(2) sector. No two- and three loop Bethe ansatz is known, except for the 
su(2) sector, where it was constructed in |jl5|. The su(2) dilatation operator was 
then extended to five loops, assuming integrability as well as BMN scaling. A five- 
loop Bethe ansatz was found experimentally, and an all-loop (asymptotic) Bethe 
ansatz was proposed |T^. For a closely related system, see |T^. Note that the 
three-loop predictions of the three-loop dilatation operator in the su(2) sector were 
recently spectacularly confirmed, for two states, by a full-fledged, rigorous field theory 



computation |T^. For a detailed review on most aspects concerning the A/" = 4 



dilatation operator, see IT^ 



One thing to note about the higher loop dilatation operators [0, |T^, is that 
they are very complicated, and appear to increase exponentially in complexity as the 
loop order increases (see e.g. appendix A.l in [0). In fact, the 5u(2|3) dilatation 



operator 0] is only known in "algorithmic" form. 

In contradistinction, the S-matrices appear to allow for a much simpler descrip- 
tion. E.g. the (asymptotic) S-matrix of [16] may be written in a very compact 
fashion (c/ ( [4.26| ) with ( |4.271) ). We shall see another example in this paper, where 
we extract the three- loop S-matrix of the su(l|l) subsector from |T^. This requires 
designing a new technique which one might term PABA: Perturbative Asymptotic 
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Bethe Ansatz. To motivate it, we derive in some, hopefully pedagogical, detail vari- 
ous one-loop Bethe ansatze at the beginning of the paper. The final result for su(l|l) 
is rather compact, while the "source", namely the two- and three-loop dilatation op- 
erator, fills at least a page. Clearly things will get worse for the operator at even 
higher loops, while we suspect the existence of a simple all-loop expression for the 
S-matrix. 

In light of the AdS/CFT correspondence it is very natural to ask about the S- 
matrix of the string sigma model. However, the sigma model believed to describe 
strings on AdS^ x is far more complicated than those theories where S-matrices 
have already been found. Still, one could try to look for similar models in order to 
derive some hints, which is an approach that has recently been followed in pO[| . 

Alternatively, one can attempt to derive useful clues from available information 
on the sigma model's spectrum. It comes from two sources. We already mentioned 

where an equation describing the classical spectrum in the su(2) sector (in gauge 
theory connotation) was derived. Recently, the procedure was also applied to the 
s[(2) sector [^. It is found that the most general classical finite-gap solution may 
be described by an algebraic curve. But does it contain information on the S-matrix, 
which is an inherently quantum concept? Here it turns out that the crucial con- 
nection does not come from the curve as such, but from the equation describing it. 
Namely, the equation may be brought into a "scattering form" which allows to draw 
conclusions about the interactions of the quantum excitations as one approaches the 
classical limit |T^,||10[. The crucial intuition for finding the correct interpretation 



comes from the fact that at one- and two-loops the classical string sigma model 
behaves very similarly ^ (but, starting from three loops, not identical [0) to 
the gauge theory when the associated operators become very "long" . On the gauge 



side one finds rather similar equations in this "thermodynamic" situation p3|, ^ 



1^ . Since the latter, which are derived from the discrete Bethe equations, do have a 
scattering interpretation, the detailed comparison allows one to also bring the string 
equations into "scattering form" . This logic gave a hint about the string's S-matrix 



in the su(2) sector |]10|, and we will apply it below, using to learn about the 
S-matrix in the s[(2) sector. 

A second important source of spectral information comes from the plane wave 
limit of strings on AdS^, x S^, where the worldsheet theory becomes free, and may be 
thus be quantized exactly [^]. Much intuition about the nature of the elementary 
excitations in this limit comes from subsequently comparing to gauge theory, i.e. from 
considering the famous BMN limit [^. Recent studies of the near-BMN limit [^ , 



28| , p9| , |30| , |31| yield information on how these excitations begin to interact. As we 
will show below, these interactions are fully consistent with the idea of an underlying 
factorized S-matrix. In fact, we will use the results of p9| , p| to derive part of 
the string S-matrix in the fermionic sector su(l|l). This is important, as it is unclear 
what may be learned from the classical sigma model for this sector. (For some recent 
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work in this direction, see |3^.) After all, fermions are inherently quantum. They 
will actually serve as a "missing link" in the present story . . . 

. . . which finishes as follows. Comparing the approximate S-matrices extracted 
from string theory in the three sectors su(2), su(l|l) and 51(2), we find a simple 
equation relating them. We then apply the relation to three-loop gauge theory, and 
write down the S-matrix in the s[(2) sector, even though we do not know the higher 
loop dilatation operator. The resulting spectrum is then successfully checked in the 
case of two-impurity states, where the spectrum is known from superconformal sym- 



metry [33|. More excitingly, we may compare to recent results of Kotikov, Lipatov, 



Onishchenko and Velizhanin |3^, who proposed three- loop exact anomalous dimen- 
sions for A/" = 4 twist-two operators. These dimensions were first found at one-loop 
from the OPE of four-point functions in , and the first few, at two loops, in [p^ 



This author does not understand the derivation of Kotikov et.al., but it is clear that 
the crucial input is a computationally highly intensive three-loop QCD field theory 
computation recently completed after a many-year effort by Moch, Vermaseren and 



Vogt |]36| (see also the comments made in [0). It involved the evaluation of (9(10 



auxiliary integrals, and the development of cutting-edge algorithmic techniques. The 
spectrum of twist-two operators derived from our S-matrix agrees with the results of 
Kotikov et.al. 

Can we do better and find the full S-matrix of the CFT and strings on AdS in 
a more direct fashion? Note that this possibility entirely hinges on whether integra- 
bility is really exact in either theory. It does not necessarily depend on whether the 
AdS/CFT duahty |^ is actually valid dynamically. We feel that the present work 



demonstrates once again that it is surely useful. 

2. One-Loop Scattering in Planar J\f = i 

2.1 The su(2) Bosonic Sector 

This sector consists of operators of the type 

Tr + . . . = Tr + . . . , (2.1) 

where J denotes an R-charge w.r.t. S0(6) and M the number of "impurities". The 
partons Z and are two out of the three complex adjoint scalars of the A/" = 4 
model. The dots indicate that we need to consider all possible orderings of them 
inside the trace, and diagonalize the set of such operators with respect to dilatation. 
This is most easily done when interpreting the dilatation operator as a spin chain 
Hamiltonian |l|]. In the spin chain interpretation L is the chain length and M the 
number of excitations. It is convenient to open up the trace and replace it by a 
quantum mechanical state on a one dimensional lattice of L sites: 

Tr (0ZZ0...0Z) ^ (2.2) 
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Let us label the sites of this lattice by a discrete coordinate x. (So in the example 
of (|2.2|) we see that we have a at a: = 1, 4, L — 1 and a Z at x = 2, 3, L). The fact 
that we originally had a trace leads to two consequences, to be distinguished. The 
first is that, since the trace links the matrix indices of the first and last parton, the 
chain has periodic boundary conditions: x = L + 1 is to be identified with x = 1. 
The second is that the cyclicity of the trace requires us to project onto states whose 
total lattice momentum is zero, see ( p^.l2 ) below. 



The planar one-loop position space Hamiltonian reads 

L L 

^0 = 5Z (1 - ^x.,.+i) = 5Z 1(1 - ■ ^cc+i) , (2.3) 

x=l x=l 

and may be expressed with the help of the permutation operator Vx,x+i which ex- 
changes the partons at sites x and x + as first noticed in It may alternatively 
be interpreted as an 5u(2) nearest neighbor spin chain if we say that Z is a spin 
up with 2-spinor (1,0) and is a spin down with 2-spinor (0,1). Then the spin 
operator at lattice site x contains the three Pauli matrices = (cr^, cr^, o"^) and the 
Hamiltonian may be written in the the form of a sum over nearest neighbor spin-spin 
interactions. 

The two-body states are defined by 

XI X2 

l*) = Ei<.,<.,<L^(^i>^2) \...ziz...z<pz...), ^^-^^ 

where xi^2 (with xi < X2) label the positions of the two particles in the background 
of the Z particles. In position space the Schrodinger equation Hq ■ l"^) = EqI"^) 
becomes 

for X2> Xi + 1 : (2.5) 
Eo ^(xi, X2) = 2 ^(xi, X2) - ^(xi - 1, X2) - ^(xi + 1, X2) + 

+2 ^(Xi, X2) - ^(Xi, X2 - 1) - ^(Xi, X2 + 1) , 

for X2 = Xi -|- 1 : (2.6) 
Eo ^(xi, X2) = 2 *(xi, X2) - ^(xi - 1, X2) - *(xi, X2 + 1) . 



This difference equation is easily solved by Bethe's ansatz for the position 
space wave function \E'(xi,X2) which reads 

^(xi,X2) = e'fi"i+'^'2"2 + S{p2,pi) e'P^"^+'P^"^ . (2.7) 

It is based on the intuition that the partons should freely evolve down the trace with 
fixed momenta pi,P2 until they hit each other at X2 = Xi-|-1. If the system is integrable 
they should then simply either pass through each other, or else exchange momenta. 
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with an amphtude given by the S-matrix S{pi,p2)- This scattering process is non- 
diffractive if the individual momenta pk are individually conserved. If we were to 
stay with the two-body problem this would of course always be true because of total 
momentum conservation. The miracle of integrability is that it remains true for the 
M-body problem. Put differently, we cannot prove integrability by only considering 
the two-body problem, but we can certainly find the S-matrix if we assume (or know) 
that the system is integrable. 

Plugging the ansatz (|2.7| ) into the Schrodinger equation (|2.5| ),( p^ one finds that 
it is indeed satisfied if, firstly, the energy is given by the dispersion law 

A/ 

= 4 sin^ (I) , (2.8) 

k=l 

with M = 2, and, secondly, if the S-matrix is given by 

gipi+ip2 _ 2e*Pi -|- 1 

In line with intuition, the dispersion law ( |2.8|) follows from the "generic" situation 
( |2.5| ) while the S-matrix is deduced from the "colliding" situation ( p.6|) . 

All this is true for arbitrary values of the momenta pk- As always in quantum 
mechanics the eigenvalues get fixed through the boundary conditions. Imposition of 
periodic boundary conditions '^{xi,X2) = '^{x2,xi + L) on the wave function (|2.71 ) 
leads to Bethe's equations: 

e^pi^ = 5(pi,p2) and e^f^^ = S(p2,Pi) • (2.10) 

Their solution leads in many cases to complex values of the momenta pk- This simply 
reflects the fact that the partons are not freely propagating, but mutually interacting 
inside a finite volume. In particular, they can form bound states. 

The principle of non-diffractive scattering, if applicable, now allows us to take a 
big leap and immediately solve the M-body problem. The total phase factor acquired 
by a parton circling around the trace should simply be given as a product of the phase 
factors due to individual collisions with all other M — 1 partons. The Bethe equations 
become 

M 

e'P^'^ = l[S{pk,p,), k = l,...,M, (2.11) 

i=i 

with the same two-body S-matrix ( p.9| )! The total energy is of course still given by 
the sum over all M local dispersion relations ( |2.8| ). 

In our gauge theory application we have to take into account the fact that the 
trace is cyclic, which means that we are only interested in the zero-momentum sector: 

M 

$^Pfc = 0. (2.12) 

k=l 
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The one-loop anomalous dimension A of the original gauge theory operators (|2.1j ) 
is then related to the energy spectrum of the zero-momentum states of the nearest- 
neighbor spin chain through 



with 



9 



87r2 



(2.13) 



where in the present case the classical scaling dimension is Aq = L. 

Finally, one commonly expresses the su(2) Bethe equations in a form familiar 
from the algebraic Bethe ansatz |^ , as in |1], ^ . Changing variables by introducing 
the so-called Bethe roots Uk = \ cot (^) , ( |2.11| ) becomes, after rewriting the S-matrix 



Uk - h 



M 

n 

3^k 



Uk — Uj + i 
Uk — Uj — i 



M. 



(2.14) 



while the momentum constraint ( |2.12| ) and the energy 

M , j M 



turn into 



n 

k=l 



Uk + 
Uk - 



and 



En 



k=l 



ui + 



(2.15) 



2.2 The su(l|l) Fermionic Sector 

This sector consists of operators of the type 



Tr + . . . = Tr V'''' ^ 



M ry.J-l 



+ ... 



(2.16) 



where J denotes an R-charge w.r.t. S0(6) and M the number of "impurities". The 
partons are for one a complex adjoint scalars Z, and secondly an adjoint gaugino if) 
(in TV = 1 connotation). In the spin chain interpretation L is the chain length and 
M the number of excitations. 

The planar one-loop Hamiltonian reads 



Hn 



Ed 

x=l 



n 



X,X + 1 J 



(2.17) 



It may be extracted from the complete one-loop J\f = 4 dilatation operator |T3[ and 
rewritten with the help of the graded permutation operator U^^x+i which exchanges 
the partons at sites x and x + 1, picking up a minus sign if the exchange involves 
two fermions ip. It was noticed in that this Hamiltonian corresponds to a free 
lattice fermion. Here we would like to put this observation into a familiar condensed 
matter context. Let us rewrite the Hamiltonian in spin chain form by expressing it 
with the help of the three Pauli matrices cx^, a^, : 



1j 



2 ^ 
X + ly 



(2.18) 



x=l 
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In spin chain language the bosonic partons Z are "spin up" spinors (1,0) and the 
fermionic partons are "spin down" spinors (0,1). We now observe that ( p.l^ ) 
is the Hamihonian of an XY spin chain in a magnetic field, which is well known 
to correspond to free lattice fermions. It is therefore even simpler than the XXX 
Heisenberg model! (The model is isotropic in the cr^, a], plane, so we could call it an 
XX spin chain.) 

The two-body states are defined by 



X\ X2 



1^) = Ei<.,<.,<L^(^i,^2) \...z^z...zijz...), ^^-^^^ 

where xi,X2 label the positions of the two gauginos in the background of the bosonic 
Z particles. In position space the Schrodinger equation Hq ■ |\&) = i?o 1^) becomes 

for X2> Xi + 1: (2.20) 

Eo ^(xi, xa) = 2 ^(xi, X2) - ^(xi - 1, Xa) - ^(xi + 1, Xa) + 

+2 ^(Xi, X2) - ^(Xi, X2 - 1) - ^(Xi, X2 + 1) , 

for X2 = Xi + 1 : (2.21) 

Eq *(Xi, X2) = 4 ^(Xi, X2) - ^(Xi - 1, X2) - ^(Xi, X2 + 1) . 

These equations are identical to the ones of the su(2) model except for the innocent 
looking replacement of a factor of 2 by a 4 when comparing ( p. 61) and (|2.21| ). 



We now make the same Bethe ansatz (|2.7|) as in the case of the su(2) sector, and 
plug it into the difference Schrodinger equations (|2.20|) , (|2.21|) . The energy eigenvalue 
(|2.8|) remains unchanged, but the result for the S-matrix is very different: 

5,u(i|i)(pi,P2) = -l. (2.22) 

So indeed the excitations behave as free fermions, and the S-matrix (|2.22|) reflects 
free fermi statistics. In fact, here the wavefunction (|2.7|) is simply a two-body Slater 
determinant. Some care has to be taken with Fermi statistics when imposing pe- 
riodic boundary conditions; we now have \E'(xi,X2) = — ^(x2,xi -|- L). The "Bethe 
equations" are then: 

e'^^^ = -5(pi,p2) = 1 and e'^^^ = -5(p2,Pi) = 1 • (2.23) 

Due to integrability we can again immediately solve the one-loop M-body problem. 
Here there simply is no scattering at all, but we need to be careful with Fermi 
statistics. The M-body equations read 

e^P^^ = l, fc = l,...,M. (2.24) 

These are of course immediately solved 



Pk = —r- , (2-25) 
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but due to fermi statistics the integer mode number Uk are required to be all distinct. 
(The wavefunction is an M x M Slater determinant and vanishes if two momenta 
coincide) . They also have to be restricted to the fundamental Brillouin zone in order 
to avoid overcounting the states. Furthermore, in gauge theory we are again only 
interested in the zero-momentum sector Yl^=i Pk = 0, which immediately translates 
into Y.kLi = 0. 

The one-loop anomalous dimension A of the original gauge theory operators 
( p.l6|) is still given by ( p.l3| ) where the classical scaling dimension is Aq = L + | M. 



2.3 The s[(2) Derivative Sector 

This sector consists of operators of the type 

Tr D^^Z^ + . . . = Tr Z'^ + . . . , (2.26) 

where J denotes an R-charge w.r.t. S0(6) and M = S the number of "impurities" 
which corresponds to one of the two spin quantum numbers^ 5* of SO (2,4). The 
partons are for one a complex adjoint scalar Z, and secondly an adjoint lightcone 
covariant derivative D. The dots indicate that we need to consider all possible 
distributions of the covariant derivatives D onto the scalars Z inside the trace, and 
diagonalize the set of such operators with respect to dilatation. In the spin chain 
interpretation L = J is the chain length and M the number of excitations. Note that, 
unlike the previous two cases, the excitations do not contribute to the length of the 
spin chain. In particular, the number of excitations may exceed the total length of 
the chain. In this context, the length is also commonly called "twist". A "twist-two" 
operator is thus a very short spin chain of length 2. 

The one-loop Hamiltonian for this sector has a global s[(2) invariance. It is non- 
polynomial, and has been derived and described in [0. In the planar limit it may 
be considered as an integrable nearest neighbor sl(2) spin — | spin chain Let us 
express the Hamiltonian H through the Hamiltonian density Ti: 

L 

Hq = ^ 'Hx^x+i ■ (2.27) 



x=l 

12. 



Denoting partons with 0,1,2 derivatives, i.e. Z, DZ, D Z by, respectively, |0), |1), 
1 2), we can read off the action of the Hamiltonian density on a neighboring pair of 
partons at lattice sites x,x + 1 from as 

n-\i,o) = |i,o) - |o, 1) n- |o, 1) = |o, i) - |i,o) 

7^-|l,l)=2|l,l)-|2,0)-|0,2) (2.28) 

H-|2,0) = ^|2,0)-|l,l)-^| 0, 2) n-\0,2) = ^ |0, 2) - 1 1, 1) - i |2, 0) , 

^We will mostly avoid using the customary letter S for the AdS spin in this paper in order 
to prevent, however unlikely, confusions with the S-matrix. Instead, we will use M (for magnon 
number). 
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where we have abbreviated Hx,x+i by H. We have not written any terms beyond 
the two-body interactions as we shall not need them for finding the S-matrix. They 
are of course important for proving the factorizability of the many-body S-matrix, 
but not for finding it once we know (or believe) that the model is integrable. This 
observation might be very useful when constructing the Bethe ansatz for this and 
other sectors at higher loops, as it should apply there as well. It also illustrates nicely 
one of the main points of this paper, namely that the S-matrix is a much simpler 
object than the dilatation operator! 
The two-body states are defined by 

\^) = Ei<x,<.,<L'^ixi,X2) \-ZiDZ)Z...Z{DZ)Z...) , 

where Xi, X2 label the positions of the two derivatives in the background of the bosonic 
Z particles. Notice the following subtle difference as compared to (p^ ) , (|2.19|) . The 



"particles" D are not occupying their own lattice sites, but sit instead on top of the 
inert particles (=holes) Z. In particular, there may be multiple-occupancy, thus now 
Xi = X2 is an allowed configuration! In condensed matter language we would say 
that the scattering is not "hard-core" . 

In position space the Schrodinger equation Hq ■ |\E') = £^0 |^) becomes 

for X2 > xi : (2.30) 
Eo ^(xi, X2) = 2 ^(xi, X2) - ^(xi - 1, X2) - ^(xi + 1, X2) + 

+2 ^(Xi, X2) - ^(Xi, X2 - 1) - ^(Xi, X2 + 1) , 

for X2 = Xi : (2-31) 
3 1 

Eo ^(Xi, X2) = - ^(Xi, X2) - ^(Xi - 1, X2) - - ^(Xi - 1, X2 - 1) + 

3 1 

+ - *(X1, X2) - ^(Xi, X2 + 1) - - ^(Xi + 1, X2 + 1) . 



Notice that the first of these equations, (2.30), is valid for both the "generic" situation 



X2 > Xi -|- 1 as well as for the nearest neighbor situation X2 = Xi -|- 1. The second 
expression describes the ultralocal on-site interaction X2 = xi of the partons. 

We make the same Bethe ansatz ( PT7| ) as in the previous two cases. Plugging it 
into the first of the two Schrodinger difference equations ( p.30|) gives as before the 



energy eigenvalue (|2.8| ), while the equation ( p.31| ), describing the collision, yields the 
S-matrix : 

gipi+jp2 _ 2e*^2 -|- 1 

Notice that it differs from the su(2) case by the exchange pi ^ P2- As before, the im- 
position of periodic boundary conditions, together with the principle of diffractionless 
scattering, leads to the Bethe equations of the form (|2.11|). 
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Changing once again variables from the momenta pk to the Bethe roots Uk 
I cot (^), (|2Tl| ) becomes 



^^a^-u^^ fc=l,...,M. (2.33) 

Uk - Uk -Uj+l 

while the momentum constraint and the energy are given as in the su(2) case by 



The one-loop anomalous dimension A of the original gauge theory operators 
( p.26|) is again given by ( |2.13|) where now the classical scaling dimension is Aq = 



L + M. 

2.4 Embedding of the Previous Sectors into su(2,2|4) 

In the previous sections we demonstrated how to find the S-matrix of the three 
simplest two-component sectors of the superspin chain of ||^ from the key principle of 
non-diffractive scattering. In the next chapter we will show, in the concrete example 
of the fermionic sector su(l|l), how to extend the method to higher loops. This 
requires the introduction of a new technique (the perturbative asymptotic Bethe 
ansatz). The method may also be applied to more than two components, but gets 
considerably more involved, as one needs to apply the so-called nested Bethe ansatz 
This is beyond the scope of the present paper. However, here we would like 
to show how, at one-loop, the above three two-component sectors are recovered as 
special cases embedded into the complete su(2,2|4) one-loop Bethe ansatz of 0. 

In the complete super spin chain we have seven types of roots. It is useful 
to visualize them with the help of a su(2,2|4) Dynkin diagram. There are several 
possible choices, and here we will consider the "Beauty" version of 0]: 

+1 

O ® O O O ® O (2.34) 

Ml M2 M3 M4 M5 Me My 

On top of the Dynkin diagram we have indicated the Dynkin labels of the repre- 
sentation corresponding to the vacuum state of the "Beauty" description (which 
corresponds to declaring the complex scalars Z to be empty lattice sites or "holes"), 
and on the bottom the number Mk of roots corresponding to the respective node. 
This way the vacuum state of a length L chain is BPS. The full equations, which can 
be found in 0, may be concisely written using the Cartan matrix corresponding to 
this Dynkin diagram. 

The su(2) M-magnon sector is of course immediately obtained by only exciting 
the central node: Mk = Sk^M. Now, for the fermionic su(l|l) chain of section is 
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not hard to verify [0] that it corresponds to the excitation pattern 

+1 

(8) O O (2.35) 

M-2 M-1 M 

Wl Vj Uk 

where we have omitted nodes that do not carry excited roots, and have introduced 
the notation m^, Vj and wi for the three types of roots. The Bethe equations of 
become 

, i \ L M , . M-1 i 

=rx^_h^lK±lrx^_hZ^!i^ i<k<M, (2.36) 

M i M-1 . M-2 i 



M-1 i 

nWl — t>,- — 7T 
— ^ 1 1 < / < M - 2 . 



Notice the absence of "self-interactions" of the fermionic roots wi. This allows us to 
eliminate^ them in the present situation with the following argument. Introduce 

M-1 , . X Af-1 , . X 

q{w):= n (^-^:. + ^) " H (^-^^■-^) " ^^.37) 

q{w) is clearly a polynomial of degree M — 2 in w. It therefore has M — 2 algebraic 
roots, which, in light of the last of the above Bethe equations, are precisely the 
fermionic roots wi. We may therefore also write 

M-2 

q{w) = i{M -1) Y[ (w-wi) . (2.38) 
1=1 

This allows us to deduce that for all 1 < j < (M — 1) 

M-2 i / i^ M-1 

TT ~Wi-2 ^ - 2) ^ TT ^'i - - ^ ^2 39) 

This however means that the second set of Bethe equations in ( |^.3(j| ) simplifies sig- 
nificantly: 

M _ _ i 

k=l "1 "A; + 2 

^The following arguments were also independently discovered by K. Zarembo (private commu- 
nication) Q. 
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What has happened is that the fermionic roots wi completely screen the "self- 
interaction" terms of the bosonic roots vj, such that these become fermionic. Clearly 
this is a general phenomenon (for this particular excitation pattern) and could be 
described pictorially as: 



(2.41) 



M-2 M-1 



M - 1 



In the case of su(l|l) we can now iterate the procedure, see diagram ( p.35|) , and 
thus derive indeed the "free" Bethe equations (|2.24| ) (as first noticed in [^) written 
in the Bethe root plane: 



Uk + 
Uk - 



k 



M. 



(2.42) 



Notice that the locations Uk of the Bethe roots are identical in the super spin chain 
and the fermionic spin chain formulation, which is far from obvious when superficially 
comparing ( ^.3(j| ) and ( ^^.42] ). In fact, this is in line with the intuition that the central 
roots Uk roots encode, via Uk = |cot the momenta of the physical excitations. 
These clearly should not depend on the description. 

Turning to the the sl(2) derivative sector of section |2.3| , we extract from 
that the excitation pattern of the states ( 2.26 ) is M4 = M, M3 = M5 = M — 1, 
M2 = Mg = M — 2 and Mi = M7 = 0. But this means that we can again apply the 
above reduction scheme as symbolized in ( |2.41| ), this time on both sides of the super 
Dynkin diagram ( p2.34 ). We thus immediately verify that 



+1 



M-2 M-1 M M-1 M-2 



O 

M 



(2.43) 



This means that now the auxiliary roots "anti-screen" the interactions of the roots 
on the central node, and we have indeed derived the Bethe equations (|2.33|) directly 
from the super chain. This result opens the interesting possibility to study the 
distribution of the auxiliary roots in the thermodynamic limit, cf |2^, as it is now 
fairly clear how to find their locations given the above discussion. 



3. Three-Loop S-Matrix for the J\f = 4 Fermionic Sector 

Now we would like to illustrate that the principle of non-diffractive scattering is 
also very powerful beyond the one-loop level. Let us apply it to the case of su(l|l) 
where the higher loop Bethe ansatz is not yet known. However, we do know the 
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Hamiltonian up to three loops in algorithmic form from the work of |14]. Rewriting 
it in spin chain form, we find the following two-loop correction^ to the XY model 
Hamiltonian ( |2.18| ) of section |2.2| : 



H, = Y, - 1) - -M'A+i - 1) + + - (3.1) 

^(^1^1 _L ^2 2 \ 3 Lrr^fn-^ J- rr"^ 

\%\^x^x+\ "T ^ x^ x+\)^ x+2 xV^ x+V^ x+2 "T ^x+V^x-^ 



x=\ 



V2} 



- + ^l+\)^l+2 - + ^l+\)^l+2 



While it certainly has not yet been rigorously proved that it corresponds to an inte- 
grable deformation of the one-loop Hamiltonian, spectral studies for small operators 
in Ijl^ are consistent with integrability in that certain tell-tale degeneracies (the so- 
called planar pairs argued to be a hallmark of integrability in p[) indeed reappear. 
Further, compelling, evidence will come from the success of the computations below, 
as they establish that the spectrum may indeed be obtained from the principle of 
factorized scattering. 

If one naively extends the approach of section |2.2| to the higher loop case one 



quickly finds that the two-body Bethe ansatz ( p.7|) for the position space wave func- 
tion \I'(xi,X2), defined in (|2.19| ), becomes inconsistent. However, we would still 
expect that the general form of the Bethe ansatz, namely a superposition of an in- 
and outgoing plane wave, is appropriate when the particles are farther apart than the 
range of the interaction, which is, in our case, the considered order of perturbation 
theory: 

^(xi, X2) ~ e'^'^"^+'P^"^ + 5(p2,Pi) e'P^"^+'^'^"^ if xi < . (3.2) 



This is Sutherland's "asymptotic" Bethe ansatz [45]. It was used by Inozemtsev to 
find the Bethe ansatz for the hyperbolic version of his spin chain [|^, and adapted 
in |jl5[ to diagonalize the three loop dilatation operator in the bosonic 5u(2) sector. 

If the ansatz ( |3.2| ) is true, the Schrodinger equation will be satisfied in the asymp- 
totic region x\ <^ X2 with the energy value given by the sum over the energies of 
the individual partons. One checks that the three-loop dispersion law of the su(l|l) 
Hamiltonian is the same as in the case of the 5u(2) chain [|T3|, and we thus have 

M 

^ = E (4 si^' (f ) - 8 9' sin^ (f ) + 32 / sin« (|) + 0(/)) , (3.3) 

k=l 



■^The three-loop piece of has been recorded (in momentum space form) in g^l; we are 
refraining from converting it into spin chain form as it is lengthy and little instructive. We found 
it simplest to use Beisert's original code, and thank him for providing it. 
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where for the time being M = 2. Furthermore, the S-matrix S{pi,p2) may be 
extracted from the asymptotics (|3.2| ), and we expect it to still be given by a pure 
phase 

^(Pi,l>2) = -e^'^^i'P^). (3.4) 

How can we find this phase factor? We will need to adapt the method of the asymp- 
totic Bethe ansatz to the present situation, since, unlike in Inozemtsev's case, here we 
currently know the Hamiltonian to three loops only. This may be done by modifying 
the "fine structure" of the wave function close to the collision point, a technique one 
might term PABA (Perturbative Asymptotic Bethe Ansatz). We make the ansatz 
(xi < X2), accurate to 0{g'^): 

- (1 + C2(pi,P2) ^2(^2-^^) + C4(P1,P2) e'P^^i+iPl^2-ie(pi,P2) 

+ 0{g') . 

Note that this form of the wave function is clearly consistent with the asymptotic 
ansatz ( p.2|) . The intuition behind (|3.5|) is that the number of powers of the coupling 
g"^ indicates the interaction range on the lattice. It should be fairly clear how to 
extend the ansatz to even higher loop order. 

Acting with the three-loop Hamiltonian^ on this wavefunction leads to difference 
equations similar to (but obviously more involved than) ( p.20| ),( p.21| ). Substituting 



the perturbative asymptotic Bethe ansatz (|3.5| ) into these equations, one finds, after 
somewhat tedious but straightforward computations, that the Schrodinger equation 
for the nearest (x2 = + 1), next-nearest {x2 = xi + 2) and next-to-next-nearest 
situation {x2 = xi + 3) may be satisfied if we carefully fine-tune the amplitude 
correction factors -82(^15^2), C'2(Pi5P2) and -84(^1,^2)5 C4(Pi,P2)- At two loops one 
finds the linear condition 

C2ip2,Pi) = i - i e'^'P'-^'P^ + e^'P'-'^'P' B2{puP2) • (3.6) 

At three loops, the functions -82(^1,^2) and 6*2(^1,^2) become fully determined, but 
their detailed form is non-universal as it depends on various gauge parameters. Fur- 
thermore, one now has a new linear constraint, similar to, but significantly more 
complicated than (|3.6|), which relates 54(^1,^2) to 6*4(^1, ^2)- We will not display it 
as it is not very instructive, and also gauge dependent. The two-body phase shift is 



^Note that that the higher loop Hamiltonians are not uniquely determined and allow for a number 
of "gauge parameters", c/ iQ. Some of these will affect the specific form of the Hamiltonian, as 
well as the wavefunction correction factors, but not the S-matrix (and therefore the spectrum). 
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then, to three-loop order, derived to be: 



d{PiiP2) = 4 sin ) sin 



Pi-P2\ . (P2\ . 

sm — + 

2 \2 ' 



(3.7) 



1 4 . fPl\ I ■ (Pl-^P2\ „ . (P1-P2 , . 

+ c/ sm ( — j I sm I 1 - 7 sm I — | + 



+ sm I I + sm 



sm I y 1 + 



There are many equivalent forms to write this phase shift; a further, interesting one 
is 

e{puP2) = (^2/ sin^ sinp2 -2/ sin^ (^) sin (2^2) + (3.8) 

+ 8g' sinpi sin^ (|) sin^ (|) + 0(/)^ - 
- (Pi ^ P2) ■ 



The final steps are again identical to the one loop case. Upon imposing periodic 
boundary conditions, taking into account Fermi statistics, the Bethe equations are 
expected to be 



M 

n 



(3.9) 



with X]fe=iPfe = 0. Now, since we already know their explicit one-loop solution from 
section we can go further and solve our equations ( |3.9| ) exactly to three-loop 
order. Taking a logarithm, we find the fundamental equation {k = 1, . . . , M) 

M 

PkL = 2Tink + ^9{pk,Pj) , (3.10) 

which is the higher loop generalization of the free fermion result ( p.25| ). Since we 
know the scattering phase shift 9{pk,Pj) exactly to three loops, cf ( |3.71 ), we may solve 
the fundamental equations recursively to 0{g^) and thus find the loop corrections 
to the free fermion momenta ( |2.25| ). Finally we substitute the obtained momenta pk 
into the expression ( p.3| ) for the energy, keeping all terms to precision 0{g^). This 
procedure then leads to the following explicit result for the, respectively, one-, two- 
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and three-loop anomalous dimensions, where -Eo,2,4 = -£'0,2,4(-^, {^«})- 

M 

^o = 4 5^sin^(^) , (3.11) 

k=l 
M 

k=l 

+ [-17) l^J (— j ^" 



A/ 



E4 = 32 5:sin«(^) + 



fc=i 



16 /7mfc\ 9 /^lA /Trfnz, — ra,) 
\ ( \ c,v.^ ( \ ( ' sm ^ ^ 1 X 



T [-17) l^J I" 



k,j=i 



L J \ L 



X (5 sin^ (^) + sin^ (^) + sin^ ^^^^AJ^ + 
16 f'^nk\ /vrnfc\ fT^ni\ /vrn^ 



fe, J, 771=1 



sm ^— X 



X sm I ^-^^ I I cos I ^ I — cos ' 



T^m ) 



Z ^ . (Tmk\ . (Tmm\ . 



k,j,m=l 



L 



. ,'2Tmj\ . f27r{nj-nk)\ , . f27i:{nj + nk) 
X sm — + sm — + sm ' 



L \ L \ L 



— 3 sm — — 3 sm ' 



L 

where the relation to the gauge theory scaling dimensions of the operators ( p.l6| ) is 
given by 

A = L + lM + g'Eo + g^E2 + g'E, + 0{g'), with / = ^- (3-12) 

We should stress that these formulas give the explicit, complete three-loop spec- 
trum of planar anomalous dimensions of all A/" = 4 operators of the form Tr Z^~'^'-'ip'^ . 
The fact that it is possible to find such a result is one of the amazing consequences 
of the higher- loop integrability of the A/" = 4 gauge theory, as first conjectured in [Q]. 
They may also be considered as an all-impurity generalization of the M = 2 formulas 
first presented in to arbitrary M . (Recall that the spectrum of "two-impurity" 
states agrees in all sectors Indeed, one checks that the expressions in ( ^.11| ) 
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reduce, after putting L = J + 1 and ni = —n2 := n, to eqs.(5.17) and (8.10) of 
One may also successfully compare to the spectrum of the first few lowest states as 
obtained by direct diagonalization of the three-loop Hamiltonian for small lengths 
L. In particular, Beisert worked out the full spectrum of highest weight states of the 
sector su(2|3) up to dimension Aq = 8.5 (see also [0, p. 143). While none of 
the su(l|l) operators are primaries in 3u(2|3), it is straightforward to find the weight 
of the corresponding primary operator, and compare its anomalous dimension with 
( p.ll|) . We present the results of this comparison in Table 1. The agreement is 
perfect. 

One may also check that the formulas ( p.ll[) explain the numerical results ob- 
tained in ||30|. Finally, it is straightforward to extract the leading 1/J correction 



to the BMN limit, for an arbitrary number of "impurities", from ( p.ll|) . The result 
agrees with near-BMN string theory at two loops, and disagrees at three, as was 
found (numerically) in |30|1. 



4. S-Matrices for Quantum Strings at Large Tension 

Let us now turn to the string side and discuss the, admittedly circumstantial, evi- 
dence that the quantum sigma model might also be described by elementary exci- 
tations living on a circle of length^ L, and whose scattering is diffract ionless. That 
is, we will try to interpret all available information on the string spectrum in the 
basic two-component sectors su(2), su(l|l) and s[(2) in the light of a (conjectured) 
factorized S-matrix. We then expect that the underlying equation describing this 
spectrum should be of the type of a fundamental equation 

M 

PkL = 2nnk + ^ 6{pk,Pj) , (4.1) 

with ^Pfc = 0, obtained by taking a logarithm on both sides of the Bethe equation. 
Here the scattering phase shift 6{pk,pj) is related to the S-matrix as 0{pk,Pj) = 
— ilog {±S{pk,Pj)) (the upper sign is for bosons and the lower for fermions), pk are 
the momenta of individual excitations, and the are quantum numbers. The total 
momentum P = Qi, the total energy E = Q2, as well as all other higher charges Qr 
are then expected to be given as linear sums over local dispersion laws qr{pk)'- 

M 

Qr = Yl <lriPk) ■ (4.2) 
k=l 

^One of the deep questions relates to the possible meaning of this "length" in the string sigma 
model. This is clearly related to the question of the nature of the elementary excitations on the 
string side. 
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Finite Length Dispersion Laws The first charge qi{pk) is the momentum p^. 



qi{pk)=Pk- (4.3) 

In we proposed for the second charge q2{Pk), i-e. for the energy per excitation, 
the following long-range dispersion law 



q2{Pk) = ^(Jl + 8g^ sin^ (f) " ^ ) ' ^^'^^ 



with 

where a/A is the string tension. This is simply a lattice version of the famous BMN 
energy formula |2^; in fact, one may show that this is the only possible lattice 
discretization of the BMN expression 0, |3^ . 

Furthermore, a study of the dispersion laws for the higher charges in the su(2) 
sector, up to five loops, suggests [|1^ that these are likely given by the expressions 

r-l 



,,(p,) = ^ — j ■ 

We expect (^4.6|) to hold for the elementary excitations in all sectors. 

BMN Limit In this limit |^ one takes L —>■ oo while keeping M = 2, 3, . . . small. 
This is a dilute gas approximation, where, except for level matching, i.e. momentum 
conservation, the excitations, in both gauge and string theory, do not feel each others 
presence. Put differently [0], we have 

e{pk,Pj)^0, (4.7) 



in ( [4.1|) , i.e. there is no scattering at all! The fundamental equation ( [4.1|) then simply 
leads to the quantization law of free non- interactive particles on a circle: pk L = 27™^. 
Substituting this solution of the "Bethe equations" into the dispersion law ( ^4.4| ) one 
then immediately recovers, via A = 1 + g"^ Q2 and using L J, the BMN formula 

A = J + X;yi + A^+^ + (^(j^) with A' = ^. (4.8) 

k=l 

where we have also already introduced the notation 5 A for the leading order 0(1/ J) 
correction to the energy, see below. 

As was first suggested in [|16], |10[, this simple picture leads to an interesting 
interpretation of the discrepancies between gauge and string theory which show up in 
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the near-BMN limit and the Frolov-Tseythn hmit [|T^, but, strikingly, apparently 
not in the strict BMN limit: The discrepancy might be due to a change in the S- 
matrix (and thus the phase shift 0{pk,Pj)) as we go from weak to strong coupling''. In 
this picture, the local dispersion laws should hold for arbitrary values of the coupling 
constant. Let us next discuss these refined limits, where the interactive but integrable 
nature of the CFT/AdS system begins to emerge. 

Near-BMN Limit Here one is interested in the first (9(1/ J) correction to the 
BMN expression ( [4.8|) . One then expects corrections^ to the free particle motion: 

What happens is that, due to finite volume effects, the dilute gas approximation 
breaks down and leading order scattering phase shifts have to be taken into account. 
If the S-matrix is known, the leading momentum shifts are immediately found from 
the fundamental equation (|47 



M 

Spk = Lj2diT^k,^nj), (4.10) 

where it is understood that we only keep terms to leading order in 1/L in the expres- 
sion 9{^nk, x^i)- Some care has to be taken as the difference between the length 
L of the system and the R-charge J begins to matter at this order. For the three 
basic two-component sectors one has L = J + uM, where, z/ = 1, |, for, respectively 
su(2), su(l|l) and s[(2). Expanding the dispersion law ( [4. 41 ) to this order, we find 
the following near-BMN energy shift 

^^-^'t^f^C-B-^MnX (4.11) 



We may then combine ( ^.lOj ) and ( [4.11| ) into the final formula 



5 A = y 5^ - eCfuk, f n,) + u {n, - n,) . (4.12) 



k,j=l 



''it will also be apparent that both limits probe the phase function and thus the S-matrix to 
leading order in small values for the individual excitation momenta. These are, unlike in finite 
length gauge theory, of order 0{1/L) in both the near-BMN and FT situation. There are indeed 
indications that the discrepancies might worsen when higher corrections (i.e. 1/J^ corrections to 
BMN and 1/J corrections to FT) are considered This is entirely consistent with our 

suspicion that the S-matrix is to be blamed. It presumably changes as one goes from weak to 
strong coupling. 

^There are some subtleties, not important for the present discussions, when some of the mode 
numbers coincide ||l|, |l^ . It was explained in detail in pl| how to deal with them, and the procedure 
immediately applies to the more general setting discussed here. 
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Expressions of this structure are indeed found in recent studies of the exact multi- 



oscillator quantization of strings in the near-plane wave geometry |2^, ^ . Below, cf 



section |4.2| , we will use this connection in reverse, and obtain important information 
on the string S-matrices from the results of 



Prolov-Tseytlin Limit In the thermodynamic limit, around the ferromagnetic 
vacuum corresponding to BPS states, the fundamental equation (|41|) is expected to 
turn into an integral equation after the rescaling p^L — > |23, |2l, Il5|, 1161 



p{(p) = 2Trn,y+ + dip' p{(p') 6{(p, if') with ip G . 



(4.13) 



This limit was originally considered in order to to compare gauge and string theory. 

Note that (^A3|) is 



following an inspiring proposal by Frolov and Tseytlin [22 



analogous to the equations (4. 9), ( 4. 10 ) of the near-BMN situation in that we expand 
the phase shift to the same leading order in 0{1/L). The only difference is that 
a large number M = 0{L) of excitations have degenerate mode numbers, which 
requires to work out how their mutual degeneracy is lifted. (This is analogous to 
the refined near-BMN situation where one also considers coinciding mode numbers.) 
In ([4.13|) p is a. convenient spectral parameter which leads to a simple form of the 
momentum p = p{p), all other charges qr{p) '■= (lr{p{p)) as well as the phase shift 
9{p,p') := 9{p{p),p{p')). The discrete "Bethe roots" pk are expected to densely 
assemble on a union of smooth contours C = Ci U C2 . . ., such that all roots on 
each component C^, carry the same quantum number n^. Note that this excludes the 
existence of this particular type of thermodynamic limit in the fermionic sector, as 
Fermi statistics does not allow coinciding quantum numbers. p{p) is the distribution 
density of the Bethe roots p, with support on C in the complex plane. It is normalized 
as 

M 

T 



— = dpp{p) 



(4.14) 



i.e. it counts the number of excitations in units of the length of the circle (however, 
c/ footnote at the beginning of this chapter). The charges ( |4.2|) are then given in the 
thermodynamic limit, after rescaling qrL"^ — > qr and QrL^ — > Qr, by 



Qr 



dpp{p)qr{p). 



(4.15) 



Equations very similar to ([4.13|) , ([4.14|) , ( |4.15[) are indeed found in recent studies of 
semi classical strings moving on the sphere or on AdS ^l|. In [|10[ we used this 



connection in reverse, and obtained important information on the string S-matrix of 
the 5u(2) sector from the results of 0. The salient points will be briefly reviewed 
in the next section 



LT| . Later we show in section |4.3| that similar results may be 
obtained for the s[{2) sector, where equations resembling Bethe equations have also 
become available [Bill. 
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4.1 The su(2) Bosonic Sector 

In an equation describing the semiclassical spectrum of strings moving with two 
large angular momenta on the five-sphere was obtained. We would expect the corre- 



sponding states to be related to gauge theory operators in the su(2) sector [^. This 
equation reads 

X + 200"^ Ex n f J f ^i.^') 



r9 



-fdx'^^ with xgC (4.16) 
Jc x-x' 



with the energy E and the momentum P given by 

E= I dx^ and P= [ dx^, (4.17) 
Jc X Jc X 

where P is quantized by an integer m: P = 2Trm, and the coupling constant uj 

appropriate for the thermodynamic limit is 

-2L^-WL^- ^^-^^^ 
(For comparison, note that cu^ is denoted by T = cu^ in [2T[].) 

The spectral equation ( 4.16|) superficially resembles the fundamental equation 
( [4.13| ) of the scattering approach, and one might be tempted to identify the spectral 
parameter x with ip in ( |4.13| ). However, the fact that an extensive quantity, namely 
the energy E, appears on the left hand side of ( 4.16 ) prevents us from interpreting 
the latter as a parton momentum p{x). A related problem is that the function cr(x) 
is found in to be normalized as 

M. = l^dxaix) (1-^) , (4-19) 

which does not allow for an interpretation of cr{x) as an excitation density as in 
( [4.14| ). However, it was shown in |T6| that both problems may be solved "in one go" 
if we change spectral parameters from x to ip according to 

^2 

= x H , with p((p) := cr(x) . (4.20) 

X 

For one, the quantity p{p) in ( [4.20| ) is then indeed normalized as in ( [4. 141) , and 
secondly, the spectral equation ( [4.16|) may now be rewritten in the form of a fun- 
damental scattering equation (|4.13| ) with the following dependence of the excitation 
momenta on the spectral parameter: 

p{p) = ^^==. (4.21) 



^In much of the hterature on the sigma model, including ||T^, the total string energy is 1 + 2lo^E. 
Here we prefer to continue to use spin chain terminology where the energy E is the eigenvalue of the 
Hamiltonian, i.e. the anomalous part of the dilatation operator divided by g^. We will also continue 
to use the same symbols in the finite L and the (rescaled) large L cases, as it should always be clear 
from the context what is meant. 
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The two-body phase shift is then found to be 



CS(^' = 2 ^o(^, + 2 5^ , (4.22) 

r=2 

where 

e^{^,^') = —^ (4.23) 
with an infinite number of additional contributions 2 9^ which may be suggestively 



written as 10 



er{^, ¥?') = u;^'- q^^^y) - qr+i{ip) Qriv')) . (4.24) 

Here the Qr are found to be given by the expressions 

<lM = ^T^T / • (4-25) 

It so turns out that these formulas are identical to the thermodynamic limit of our 
lattice dispersion laws (|4.6|) when expressed, via ( [4.21|) , as functions of the spectral 
parameter ip. This is fascinating, as they simply follow from rewriting the semiclas- 
sical equations of |^ in a form that unveils the underlying scattering processes! 

What is more, recovering the thermodynamic fundamental scattering equation 
from classical string theory suggests a very natural way to rediscretize it and guess 
the large tension, small momentum (i.e. large length L) S-matrix of the quantum 



sigma model |I0|. Here we recall that the proposed long-range S-matrix of the 
weak-coupling gauge theory in the su(2) sector is 

Its logarithm yields, in the thermodynamic limit, the leading contribution Oq, see 
( [4.23|) , to the string theory two-body scattering kernel. The all-order expression for 
the lattice phase function (p{pk) is conjectured |jl6[| to be 



V'(j>.) = |cot(|)^l + 8^2sin^(|), (4.27) 

and yields in the thermodynamic limit, upon inversion, (|4.21 ). It is then very natural 
to include the additional scattering terms 9r of ( [4. 241) as a "dressing factor" to the 
"bare" S-matrix ([4.26|) and write down an ansatz for the string S-matrix [p!0|] : 



^u(2)(P^,P.)-^(^^)_^(^^.)_,lle , (4.28 
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while replacing (1), via ( 4.18| ), the continuum coupling lJ^ by the lattice coupling 



and (2) the continuum dispersion laws (|4.25|) by the lattice laws 



OriPk,Pj) = (^^■(P'^) 1r+liPj) - (lr+l{Pk) qriPj)) ■ (4.29) 

One now checks, with the help of (|4.12|) , agreement with the near-BMN multi- 



impurity results of |29, 31|, as was first done in |10 



4.2 The su(l|l) Fermionic Sector 

In the fermionic sector no "ferromagnetic" thermodynamic limit similar to the one 
discussed in the last section exists. Due to Fermi statistics the mode numbers of all 
excitations are distinct. This prevents the Bethe roots from condensing onto smooth 
cuts in the complex spectral parameter plane, where the mode number has to stay 
constant along each contour. Luckily there is an alternative way to deduce informa- 
tion on the string S-matrix of this sector. In a beautiful paper Callan, McLoughlin 



and Swanson |]2D| studied the spectrum of three elementary string excitations in the 
near-BMN limit. It so turns out that this is precisely what we need! In section ^ we 
argued that, given integrability, it suffices to carefully solve the two-body problem in 
order to deduce the S-matrix. It is however important that the two excitations are 
"off-shell", i.e. they need to be capable of carrying arbitrary momenta. Now, since 
the level matching condition in string theory enforces the total momentum conserva- 



tion of all excitations, the solution of the two-impurity problem |2^ does not contain 
enough information as the absolute values of the two associated momenta are nec- 
essarily equal. This is no longer the case if we decompose the three-body scattering 
into a series of two-body processes. And indeed, the principle of factorized scatter- 
ing then immediately yields the solution of the M-body problem, as was intuitively 
and correctly understood in a nice follow-up paper by McLoughlin and Swanson 
Rewritten in our present notations, they find (c/eq.(3.30) in [^) the following 



energy shift for M excitations in the su(l|l) sector: 

SA = -- >^ I / ' -2n,n, \ . (4.30) 

/! + \'nlJl + y 



This double sum may be rearranged as 

M 

6A = X' y , X (4.31) 

k% V^TX^k 



X I 



Y^l + A'n| -l\-lnk {^^l + X'n] - l) + \ (n, - n^) 
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By comparing to our general formula ( 4.12|) (here u = j as appropriate for the 



fermionic sector), we may then read of the phase shift, accurate to leading order in 
small pk, Pj, as 



-0i{Pk,Pj] 



9 

-^{Pkq2{Pj) -q2{Pk)Pj) 



(4.32) 



where we have used ( [4.4|) (for small pk) and employed the notation 6r{pk,Pj), intro- 
duced in ( [4.29|) , with r = 1. This then leads to the following simple, approximate 
S-matrix encoding the near-BMN string physics in the fermi sector: 



_p-ifi(pfe,pj) 



(4.33) 



4.3 The 5[(2) Derivative Sector 

In the case of the third basic two-component sector si{2) we have a choice for extract- 
ing the S-matrix from known string theory results. We could either proceed as in the 
last section using again (|4.12|) in conjunction with the multi-impurity near-BMN 
results of [^, |31| . Alternatively we can apply once more the logic of section |4.1| and 
test them on the semiclassical string sigma model results of Kazakov and Zarembo 
which have recently become available Interestingly, both procedures lead to 

the same final result, showing once more the close connection between the leading 
curvature corrections to quantum strings in a plane wave geometry, and semiclassical 
strings on the curved geometry AdS^ x S^. Let us proceed in the second fashion. 



which will serve as a nice check on the ideas presented in section 



For semiclassical strings rotating with one large angular momentum on the five- 



sphere and one large spin on AdS, the bootstrap equation derived in from the 
monodromy matrix of the classically integrable string sigma model reads 



X 



2 7rn„ 



dx 



with 



(4.34) 



where the function cr(x) is normalized exactly as in the case according to formula 
( [4.19|) (now M is the spin quantum number on AdSs, whereas previously M was one 
of the two angular momenta on S^, see also ||2^). The equation ( [4.34|) is manifestly 
very similar to the corresponding equation ([4.16|) of the su(2) sector. Again an 
extensive quantity appears on the left side of this equation, preventing its naive 
interpretation as a fundamental scattering equation: This time it is not the energy 
E but the total momentum P, see ( [4.171 ). However, since the normalization conditions 
( [4. 191 ) are identical in the two cases, it is natural to once more apply the same change 
of spectral parameter ( [4.20|) . This converts (|4.34| ) to the form of a fundamental 
scattering equation, see ([4.13|) 



p{lp) = 2 71711, - f dip' p{ip') 



- + 2^^,(y. 




(4.35) 
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where G Cj/, the momentum is parametrized as before through (|4.21|) , the density 
p{ip) is consistently normalized as in (^4.14|) , and the phase shifts 6r{(p,(p') are given 
in ( |4.24| ). We may then read off the continuum two-body scattering phase shift of 
the s[(2) sector as 

oo 

C(2T(V'' = -2 ^o(v^, ^riv^ ' (4-36) 

r=l 

which differs from the one of the su(2) sector ( [4.22| ) by (1) an overall minus sign and 
(2) the fact that the sum starts at r = 1 instead of r = 2. 

As a consistency check we may verify that the near-BMN physics is properly 
reproduced. Replacing the shift (|4.36| ) by its discrete, small pk (i.e. large L) version 
in the, by now, familiar fashion, we find 

2 °° 

OiPk.Pj) ^ ^ - 2 V e^ipk,?,) , (4.37) 



where the lattice expressions ( [4.27|) and (|4.29| ) should be used (replacing momenta 



by their leading 0{1/L) approximations). If we now insert (|4.37|) into the formula 
( [4.12|) for the energy shift (with z/ = for the s[(2) sector) we reproduce, after a 



short calculation, indeed the sl(2) multi-impurity result of [^, ^ 

Uk rij — X' nl n'j 




+ Uk Uj . (4.38) 



c/ equation (3.24) in [^]. The case of coinciding mode numbers is dealt with exactly 



as m 



1T0[] , and immediately reproduces equation (3.23) of |]3l| . 
Finally we may write the (large tension, small momentum) S-matrix for the s[(2) 
sector suggested by the phase shift ([4.37|) : 



^5? p.) - ^i^i fl e-"^^^-'^^^ • (4.39) 

Thus the approximate string S-matrix in the si{2) sector differs from the one of the 
su(2) sector (|4.28| ) of [T^ by (1) complex conjugating the explicit factors of 2 — > — z 



and (2) by one additional factor with r = 1. 

5. Three-Loop S-Matrix for the = 4 Derivative Sector 

Observe the following striking relationship between the conjectured large string 
tension, small momentum S-matrices ( 4.28| ),( ^7S^ ),( ^7SU ) of the three basic two- 



component sectors: 

Ssl{2) = 5'su(i|i) 5'^u(2) 'S'su(lll) • (5.1) 



-27- 



The simplicity of the relation suggests a purely group-theoretical explanation. 
In fact, it appears to result from an "inversion" of the Dynkin diagram (p.34|) where 
we now place the representation —1 onto the central node; i.e. we exchange the AdS 
and sectors. 

Since the relationship seems to be based on symmetry alone it is very 

reasonable to expect that it also holds at weak coupling, i.e. in the gauge theory. We 
are thus led to the following Bethe ansatz for the weak-coupling s[(2) sector at finite 
L = J and finite M = S (with A; = 1, . . . , M): 



M 



e 



ipkL 



where (p{p) is given in ( [4. 271) and 0{pk,Pj) is currently known to three-loop precision, 
cf (P) or (pD. 

A first, non-trivial test of this ansatz may be performed immediately. We know 
that all two-impurity states must agree for arbitrary, finite R-charge J ||33|. Let us 



compare the M = 2 Bethe ansatze for the su(2) sectors O, fB] and the su(l|l 



sectors ( cf section |2l^ ) : 



su(2) : e*f(-^+2) = ^M±i 5u(l|l) : e'^^'+'^ = e^^^^'-^) , (5.3) 

from which we conclude that we must have to all orders in the coupling constant g 

gie(p -p) ^ g-ip '^(P) + t _ (^5 4^ 

The reason is that we have the same dispersion relation in all sectors, and the mo- 
menta pi '■= p, P2 '■= —p for two impurities must therefore agree in different sectors. 
We easily check that the three-loop phase shift obtained in section |2.2| satisfies ( |5.4| ) 



to 0{g ). One now immediately verifies that the ansatz (|5.2|) for the sl{2) sector 



where the length is L = J, results in a Bethe equation for p which is entirely equiv- 
alent to (|0| ): 

sl(2) : e'P' = e^inv,-P) (5.5) 

vi.p) + I 



Now, naively we would expect the Bethe ansatz (|5.2|) to be asymptotic only, that 
is we may a priori not expect that it properly diagonalizes operators which are shorter 
than the range of interaction. These are precisely the twist-two operators at two and 
three loops, and the twist three operators at three loops. However, the consistency 
check just presented actually yields a further, crucial hint. We notice that our si{2) 
ansatz also works at least up to three loops for the M = 2 states with J = 2 and J = 
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3, cf (|5.5| ). But these are, respectively, length two and length three operators. This 
observation leads to the expectation that in fact all sl{2) operators are diagonalized 
by the Bethe ansatz ( [5 .21 )! This allows us to test our ansatz since, luckily, the 
three-loop anomalous dimensions of twist-two operators in the Af = 4 gauge theory 
have recently become available. Kotikov, Lipatov, Onishchenko and Velizhanin p5| 



were able to extract them, under some unproven but astute assumptions, from an 
impressive, rigorous computation of three-loop anomalous dimensions in QCD by 
Moch, Vermaseren and Vogt []36 |. 



We have not worked out the explicit three-loop spectrum of twist two operators 
from (|5.2| ) for arbitrary even AdS spin M (in the sl{2) sector S = M has to be even 
for twist-two operators), but certainly expect that this will reproduce the general 



formulas derived in ||35[ (c/ equations (10), (11), (12) of that paper, where the notation 
J = M is used), which read: 

Eo(M)=4Si, (5.6) 

E^iM) = -4 (S3 + S_3 - 2 S_2,i + 2 Si (S2 + S_2)) , 

E4{M) = -8 (2 S_3 S2 - S5 - 2 S_2 S3 - 3 S_5 + 24 S_2,i,i,i 

+ 6 (S_4^i + S_3^2 + S-2,3) — 12 (S_3^i^i + S_2,l,2 + S_2,2,l) 

- (S2 + 2 S?) (3 S_3 + S3 - 2 S_2,i) - Si (8 S_4 + ^-2 
+ 4 S2 S_2 + 2 S^ + 3 S4 - 12 S_3,i - 10 S_2,2 + 16 S_2,i,i)) , 
and the harmonic sums T,±a,b,c,--- '■= ^±a,b,c,---{^) defined recursively {a,b,c > 0) 

S±.(M) = V^^, S±,,,,,,..(M) = V^^S,,,,..(m). 

m=l m=l 



We checked explicitly that our ansatz (|5^) reproduces the result predicted in ([57 
in the cases M = 2,4, 6, 8; see also Table 2. 

We would however like to stress that the ansatz ( |5.2| ) is expected to also properly 
reproduce the three-loop anomalous dimensions for operators of arbitrary twist and 
spin. As far as we know no non-trivial twist-three two-loop, let alone three-loop, 
anomalous dimensions seem to be known in A/" = 4 gauge theory from a field theory 
computation to date. For example, for the simplest twist-three field not in the two- 
impurity supermultiplet, namely Ti D^Z^ + . . ., our ansatz predicts a paired state 
with three-loop energy 

„ 15 225 o 3195 a , a. , , 

It would be exciting if this prediction could be checked by a traditional Feynman 
diagram computation^^ 

After the completion of this manuscript we were informed by B. Eden that the two-loop part 



of the prediction (5.7) may indeed be confirmed by a full-fledged field theory computation W9l. 
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Finally we may answer a question raised in the paper [^] about whether the 
discrepancies between string and gauge theory might already appear in the s[(2) 
sector at two instead of at three loops. Comparing the S-matrix entering the r.h.s. of 



( ^■2|) (c/also (|6.4|) ) and the string S-matrix ( 4.39 ) by expanding in small momenta 
Pk, Pj of order 0{1/L) one finds, 

-t (log Sl^i^ip,, pj) - log p,)) ^^pI {pk - Pj)p] + 0{g') . (5.8) 

This proves, incidentally for both the near-BMN limit as well as the Frolov-Tseytlin 
limit, two-loop agreement and three-loop disagreement in the 5l(2) sector, in full 
analogy with the su(2) and su(l|l) sectors. 

6. Summary and Musing 

Let us collect once more the long-range, asymptotic Bethe ansatze we proposed for 
the basic two-component sectors su(2), su(l|l) and s[(2) of A/" = 4 gauge theory. 
They are all written in the factorized scattering form 

M 

e^P.i^ -Q ±s{pk,pj), k = l,...,M, (6.1) 

where the upper sign is for bosons and the lower for fermions. The three S-matrices 
are 

= (6.2) 

SI:;IT\iMp^) = -e''^''''^K (6.3) 
<^(Pfc) - (pjPj) - i 
f{Pk) - <f{Pj) + i 



s^^^{Pk.p,) = ^ ~ 7 ] e2^^(^-^^) . (6.4) 



The phase function ^p{pk) is currently known to five loops, and its all-loop conjecture 
is (|4.27| ) The phase 0{pk,Pj) was worked out in section | to three-loop order, cf 
( p.7|) , (|3.8|) . It would be fascinating to find its all-loop form. 

On the string side, it was found, for su(2) in , and in section ^ for the other two 
sectors, that the near-BMN and Frolov-Tseytlin physics is, for large L (i.e. momenta 
of order C(l/L)) reproduced to all orders in by the S-matrices 

st^ij^k, p,) - ^;^^|-^;^^|+^ n ^"'^^"^^^ > (6-5) 

C;i|i)(P'^'P.)^-e-^'^^^'='^^7 (6.6) 



r=l 
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Finally, it was demonstrated in that to at least five-loop order a spin chain exists 



whose S-matrix Ssu{2) is exactly given by the right hand side of ( |6.5|) . Interestingly this 
model is part of the general three- loop su(2|3) spin chain of which also contains 
the su(l|l) subsector. Therefore supersymmetry predicts that also a fermionic spin 
chain should exist which reproduces the string and not the gauge theory scattering. 
However, as opposed to (|6.5D, the S-matrix ( |S.6D can unfortunately not, as it stands, 
correspond to the one of a quantum spin chain. The reason is that the shift function 
Or{pk,Pj) is not periodic in the momenta Pk,Pj for r = 1. Excitingly, one may 
nevertheless find evidence that the same deformation which turns the gauge theory 
spin chain into the string theory spin chain in the su(2) sector yields for the su(l|l) 
sector an S-matrix S'su(i|i) which is, to three loop order, "dressed" by the same 
multiplicative factor exp2i92{pk,Pj) as in the su(2) case. In fact, we have checked, 
by including the corresponding deformation parameter (denoted by C4 in |^0[, where 
C4 is the gauge case and C4 = 1 the string case) into the calculations of section |^ that 
the deformation leads to the following additional phase shift for su(l|l) 

Sd{pi,P2) = -8c4 5f^sin^ (^y) (sin Pi - smp2^ sin^ (yj , (6.8) 

to be added to (|3.8| ). (Note that in the string case C4 = 1 it exactly cancels the last 
term in ( |3.8|) .) This is possible since the deformation parameter appears in the full 
three- loop su(2|3) vertex of |]14|. It is therefore reasonable to speculate that a spin 
chain might exists which reproduces the string results for all sectors and, maybe, all 
of psu(2, 2|4). For the two-component sectors discussed in this paper this spin chain 
might have an S-matrix which is related to the weak-coupling, asymptotic gauge 
theory S-matrix by the same dressing factor: 

SZl^)iPk.Pj) = SZJ^iPk.Pj) Sipk^Pj) , (6.9) 

s:zh^Pk,P,) = si::;T\i)iPk,p,)siPk,p,), 
s:;S)iPk,pj) = s'^^{pk.p,)s{pk^p,). 

where the dressing factor would be approximately given by the product 

00 

^(pfc,p,)~ JJe^^^'-^^-P^). (6.10) 

r=2 

Hopefully such a dressing factor will appear in the gauge theory when we go from 
weak to strong coupling, e.g. through the "wrapping effects" discussed in |l^. Or do 
we have to replace the AdS/CFT "correspondence" by an AdS/CFT "similarity"? 



^^It should be clear from the perturbative asymptotic Bethe ansatz mtroduced m this paper that 
the latter will always yield a periodic S-matrix. The same caveat therefore applies to the s[(2) 



string S-matrix (5.7). 
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It would be exceedingly important to test these ideas in refined situations such 
as corrections to the near-BMN limit, and 1/L corrections to spinning strings. 
For first steps in these directions see, respectively, |^ and [0, A dressing factor 
as in (|6.9|) should then also explain the famous strong coupling behavior. First 
evidence that this might happen was found in |]T0[. It would also be very interesting 
to reproduce other semiclassical limits, such as the large spin limit of [B2|, from the 
S-matrices. 
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Table 1: The three-loop spectrum of the first few excited states of the fermionic subsector 
5u(l|l) as computed from the explicit solution ( 3.11| ). The latter is derived from the pertur- 
bative Bethe ansatz. The mode numbers entering the fundamental equation are indicated. 
The eigenvalues are either rational or algebraic; in the latter case we have given them to 
five digit precision. In order to compare with the results of direct matrix diagonalization 
ll^ (see also p. 143 of [^]) we have indicated the labels of the corresponding highest 



m 



weight state of su(2|3) (no su(l|l) state is primary in su(2|3)). The agreement is perfect, 
and establishes perturbative factorized scattering in the fermionic sector. Note that the 
agreement persists up to large magnon densities, e.g. M = 4 magnons in a length L = 5 
spin chain. P denotes the parity of a state. 
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Table 2: The three-loop spectrum of the first few highest weight states of the derivative 
subsector s[(2) as derived from the perturbative Bethe ansatz. The mode numbers entering 
the fundamental equation are indicated. The eigenvalues are either rational or algebraic; 
in the latter case we have given them to five digit precision. The one-loop spectrum was 
found by direct matrix diagonalization in ||T^ (see also p. 84 of |jl^). At the time of writing 
this spectrum cannot be checked against direct diagonalization since the 5[(2) dilatation 
operator is not yet known beyond one loop. The bold-faced values with L = 2 correspond 
to twist-two operators and agree perfectly with results of Kotikov, Lipatov, Onishchenko 
and Velizhanin. The values of two-magnon states M = 2 were previously known from 
superconformal invariance. The three-loop results for the cases L = 2, 3 and M = 2 were 
first predicted in [Q] and confirmed in a field theory computation in [18|. The two- and 
three loop anomalous dimensions of all twist-L operators with L > 2 and M > 2 are 
new predictions. The parity of a state is denoted by P. The case M = 4, L = 4 with 
mode numbers ±(1, 1, 1, 1) is interesting as it is the simplest example of a paired state in 
s[(2) where all Bethe roots are either strictly positive or strictly negative real. Such states 



with "winding number" were considered, in the thermodynamic limit, in p3| , 21], and it is 
reassuring to find them in the finite spin chain. 
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